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Abstract 

We propose a new class of dark matter models with unusual phenomenology. What is ordinary 
about our models is that dark matter particles are WIMPs, they are weakly coupled to the Standard 
Model and have weak scale masses. What is unusual is that they come in multiplets of a new “dark” 
non-Abelian gauge group with milli-weak coupling. The massless dark gluons of this dark gauge 
group contribute to the energy density of the universe as a form of weakly self-interacting dark 
radiation. In this paper we explore the consequences of having i.) dark matter in multiplets ii.) 
self-interacting dark radiation and iii.) dark matter which is weakly coupled to dark radiation. We 
find that i.) dark matter cross sections are modified by multiplicity factors which have signihcant 
consequences for collider searches and indirect detection, ii.) dark gluons have thermal abundances 
which affect the CMB as dark radiation. Unlike additional massless neutrino species the dark 
gluons are interacting and have vanishing viscosity and iii.) the coupling of dark radiation to 
dark matter represents a new mechanism for damping the large scale structure power spectrum. 
A combination of additional radiation and slightly damped structure is interesting because it can 
remove tensions between global ACDM fits from the CMB and direct measurements of the Hubble 
expansion rate (Hq) and large scale structure (erg). 
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I. INTRODUCTION 


Dark matter (DM) makes up more than 80% of the matter content of our universe [T]. 
Despite the overwhelming quantitative evidence for cold dark matter (CDM) from its grav¬ 
itational interactions, and its inclusion in the standard ACDM cosmological model, we only 
have very limited information about possible non-gravitational interactions of DM. Existing 
upper bounds on interactions between the DM and the Standard Model (SM) and self¬ 
interactions of the DM leave a broad range of possibilities for models of DM. Viable models 
may include a whole dark sector consisting of different species of particles which interact 
only with themselves and with DM (see e.g. pm). Furthermore, there is some tension be¬ 
tween the best £t parameters in ACDM as inferred from the cosmic microwave background 
(CMB) [U, and direct measurements of Hq using standard candles and ag using galaxy clus¬ 
ters [3H7]). While it is premature to draw dehnite conclusions about these discrepancies, 
they motivate the study of interacting dark sectors which go beyond the ordinary cold dark 
matter of ACDM. 

In this work, we propose a novel scenario for an interacting dark sector which predicts 
several important modihcations to standard CDM phenomenology. Our dark sector includes 
a weakly interacting massive particle (WIMP) as the dark matter. As usual, weak interac¬ 
tions between the DM and the SM keep the two sectors in equilibrium at high temperatures 
and set the abundance of DM through thermal freeze-out. The non-trivial new ingredient 
of our scenario is that the DM particles are charged under a new non-Abelian “dark” gauge 
group. ^ Thus our dark sector consists of a degenerate multiplet of massive DM particles 
coupled to “dark radiation” in the form of massless non-Abelian gauge bosons. Observations 
of large scale structure tightly constrain long-range interactions of DM with dark radiation 
and require very small gauge couplings^, Qd < 10“^, so that the “dark” conhnement temper¬ 
ature is well below the current CMB temperature.^ Thus conhnement in the dark sector is 

^ New visible particles charged under a new non-Abelian gauge group with macroscopic sized confinement 
scale have been proposed in the past under the name of theta-particles [HHin] and more recently as 
quirks m- 

^ Such couplings are too small to have appreciable effects on the Bullet Cluster, the shape of galactic dark 
matter halos or DM protohalos. They are also too small to have any bearing on the “core vs cusp” or 
“missing satellites” problems. 

^ For recent discussions of non-Abelian dark sectors in a different context see [IMI]. 
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irrelevant, and the “dark gluons” are a thermal bath of weakly coupled massless particles 
which are well described as a perfect fluid. Note that this is distinct from neutrinos which 
free-stream and are therefore not a perfect fluid. The two types of dark radiation can be 
distinguished through their imprints on the CMB UHlEn]. 

The general features of our dark sector lead to several distinct observable effects which 
can be broadly classihed as due to (h) multiplicity of the dark matter, [ii.) dark radiation 
in the form of the dark gluon fluid, and {in.) interactions between DM and dark radiation. 
We briefly summarize the salient features of each of the three effects here. Details will be 
provided in subsequent sections. 

• DM multiplicity: Being charged under an unbroken non-Abelian gauge group, DM 
particles come in degenerate multiplets of N dark colors. This multiplicity leads to 
several easy-to-understand but important differences from conventional WIMPs. First 
off, DM annihilation cross sections (into SM particles) are suppressed by 1/A^ from 
averaging over initial states. Assuming thermal freeze-out, we predict that our dark 
matter mass is smaller by a factor of l/y/N than in models without dark matter 
multiplicity. Similarly, indirect detection bounds from dark matter annihilation into 
photons near the galactic center are relaxed because of the 1/A^ in the annihilation 
cross section. Second, DM pair production cross sections at colliders are enhanced 
by a factor of N from the multiplicity of possible hnal states. Third, direct detection 
bounds for DM are unchanged because the scattering cross section for DM particles 
off nuclei does not contain multiplicity factors. We discuss these effects, which would 
also follow from a non-Abelian global symmetry, in Section 3. 

• Dark radiation: Through their interactions with the DM, dark gluons come into ther¬ 
mal equilibrium with the SM in the early universe. After freeze-out of the DM the 
dark gluons decouple from the SM. They maintain a thermal distribution but end up 
cooling relative to the CMB because the photons are heated by absorbing the entropy 
contained in massive SM particles. The dark gluons are observable through their con¬ 
tribution to the radiation density in the universe, an effect which is conventionally 
expressed in terms of an effective number of neutrino species (see 121 for a recent 
discussion of bounds on new relativistic degrees of freedom from the CMB). We hnd 
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ANeff 0.07(iV^ — 1) for the iV^ — 1 gluons of an SU{N)ci gauge group. Bounds on 
Neff from Planck and from nucleosynthesis give iV < 4 at 95% conhdence level. An 
increase in the radiation density also shifts the best £t values of other cosmological pa¬ 
rameters from the CMB. The two most significant shifts are increases in the predicted 
values of Ho and as- An increase of Ho would remove the tension between CMB hts 
and direct measurements of Ho- An increase in as worsens the tension between CMB 
data and direct observations of as from large scale structure. However, as discussed 
in the next item, interactions between DM and dark gluons can remove this tension. 
Finally, the fact that dark gluons are self-interacting allows one to distinguish them 
from additional neutrino species. Dark gluons are well described by a perfect fluid 
with no viscosity and therefore lead to less damping of the CMB power spectrum in 
comparison to extra free-streaming relativistic fluids. Future CMB data will be precise 
enough to determine the viscosity of any signihcant new component of dark radiation 
and can therefore distinguish between additional neutrinos and dark gluons [TH]. We 
elaborate on the phenomenology of our dark radiation in Section 4. 

• DM - dark radiation interactions: Scattering of dark radiation (DR) off DM introduces 
a drag force between the non-relativistic DM fluid and the relativistic radiation. This 
drag suppresses gravitational clustering and is therefore observable in the matter power 
spectrum. Allowing for matter perturbations to grow at least approximately as in 
ACDM gives a conservative upper bound of ad ^ 10“®. What we hnd particularly 
interesting about our scenario is that the momentum transfer cross section for DM - 
DR scales with temperature like the Hubble parameter during radiation domination. 
This means that the effect of the radiation remains equally important throughout 
radiation domination, leading to a smooth reduction of the power spectrum. This is 
in contrast with the more frequently studied case where the DM interactions freeze 
out at a critical temperature, leading to a sharp cutoff in the power spectrum at 
small scales [H I2I1 - I23] . We show that a very small and smooth reduction of the power 
spectrum may resolve the tension in the indirect determination of as from Planck data 
and direct measurements of the power spectrum at large scales (see e.g. EEIE]).' 

^ We thank J. Lesgourgues for alerting us to the discrepancy between determinations of ug from fits to the 
CMB and from direct measurements. 
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We discuss DM-DR interactions in Section 5. 


Most of the features associated with non-Abelian dark matter are expected to hold in 
a broad class of models irrespective of the specihc interactions between the DM and the 
SM. However, in the interest of concreteness and for clarity of presentation we focus on a 
specihc example model which is described in Section 2. The DM in the model is a “wino- 
like” SU{2)weak triplet Dirac fermion which transforms as a fundamental under a dark 
SU{N)d gauge group. 


II. THE MODEL AND DM THERMAL RELIC ABUNDANCE 


For concreteness we focus on a specihc realization of a non-Abelian dark matter. We 
take the dark matter particle to be a Dirac fermion in the (1, 3)o representation of the 
Standard Model (S't/(3), S't/(2))f/(i) gauge group and in the fundamental representation of 
dark SU{N)d . It has an SU{2)weak preserving Dirac mass and no additional interac¬ 
tions with the Standard Model besides gauge interactions. Electroweak symmetry breaking 
leads to a loop-induced mass splitting between the charged and neutral components of the 
SU{2)y^eak triplet. In the limit S> M\y the splitting is independent of and given 
by [SUES] 

6M^ = M^± - M^o ^ 0.16 GeV. (1) 

The neutral component is the lightest particle which carries fundamental SU{N)d charge 
and is therefore stable, it is our DM candidate. From the perspective of SM interactions the 
dark SU (A)d is simply a global symmetry, so that our DM acts like N identical copies of a 
weakly interacting particle with Dirac mass M^. In the following, we determine by 
requiring that the correct DM abundance is obtained from thermal freeze-out. 

At the time of DM chemical freeze-out T ~ M^/26 S> 6M^ so that we can ignore the 
mass splitting between the different components of when calculating relic abundances. 
Once the temperature drops below SM^, the charged decay to x° plus standard model 
particles. Thus the (co-moving) number density of x° particles at low temperatures is simply 
given by summing over all components of x“ at DM freeze-out, i.e. This 

explains why we can use the abundance calculation for the whole triplet x“ in order to hnd 
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the DM relic abundance. 


We will show in subsequent sections that Cosmology requires the dark gauge coupling to 
be much smaller than the SM gauge couplings. Therefore the relevant interactions for the 
thermal relic calculations are of the form yy —)■ SM SM and independent of the dark gauge 
coupling. Figure [T] shows some of the relevant diagrams for DM annihilating to SM particles. 
In the limit the thermally averaged effective annihilation cross-section is 

/ ^ 1 374 

~ m 1927rM2 ’ ^ ^ 

X 

where g 2 is the SU{2)^eak gauge coupling. This cross-section differs from the standard 
SU{2)^eak triplet “wino” [2S] by the extra 1/2N factor, which comes from the multiplic¬ 
ity associated with the SU{N)d representation and the fact that y is a 4-component Dirac 
fermion instead of a 2-component Majorana fermion. The l/2iV factor can be easily under¬ 
stood from the fact that any given DM particle carries an dark color charge and can only 
annihilate if it hnds an anti-particle with the corresponding anti-dark color, thus reducing 
the color averaged annihilation cross-section. 


vwv"' 

ww"' 
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FIG. 1: Annihilation of dark matter into SU{2)u,eak gauge bosons or SM fermions. Since the 
mass splitting between members of SU{2)y^^ak multiplets is small compared with energy trans¬ 
fer in the annihilation diagrams (i.e. twice the y mass) we compute the co-annihilation of full 
SU{2)u,eak multiplets and ignore the mass splittings. 

If the DM abundance is set by thermal freeze-out, then - to a good approximation - the DM 
mass density today depends on the DM mass only through its annihilation cross-section [27] . 
Therefore, holding the cross section in Eq. [^hxed at the Cosmologically preferred value we 
see that the mass required to get the correct relic abundance decreases as the square root 
of N. In Table [I] we give the mass of DM for different values of N using the tree level 
annihilation cross-section. The masses are signihcantly lower than for the usual “wino” case 
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where the tree level formula predicts a mass of 2.4 TeV. ® 


N = 2 

N = 3 

V = 4 

N = 5 

Generic N 

1.2 TeV 

1.0 TeV 

0.9 TeV 

0.8 TeV 

~ 2A/V^ TeV 


TABLE I: Dark Matter masses required to get the correct thermal abundance as a function of N. 


III. DARK MATTER MULTIPLICITY AND EXPERIMENTAL SEARCHES 


As will be discussed in subsequent Sections, astrophysical and cosmological constraints 
require very small dark gauge couplings of order ad < 10“®. This is too small to have 
observable effects for production of DM at colliders, or for direct or indirect detection of 
DM. However, it is important to take into account that DM particles are in multiplets of 
the dark SU{N)d , i.e. there are multiplicity factors associated with DM processes. In the 
previous Section we already discussed that the predicted mass for our thermal relic DM is 
reduced by a multiplicity factor of l/\/2N, which can lead to large changes in sensitivity 
of experimental searches. Similar multiplicity/color factors also appear in cross-sections 
associated with DM detection. Figure [^shows the color factors associated with the different 
kinds of DM searches. In what follows we briefly describe the effects of SU{N)d multiplicity 
for direct and indirect detection and also for collider searches for DM: 


A. Direct detection 


For direct detection there is no color factor associated with the multiplicity, thus the only 
change comes from dark matter being lighter. The spin-independent cross-section for dark 
matter scattering of the nucleus is approximately 10“^^ cm^, and independent of the DM 
mass as long as 3> Mw [22] • This cross-section is an order of magnitude smaller than the 
projected sensitivity of the next generation direct detection experiments However, in 

® In the case of the “wino”, the preferred mass is close to a Sommerfeld resonance of the weak interactions 
and one must account for Sommerfeld enhancement to obtain M^i^o — 3 TeV |2H]- In our case, the DM 
mass is safely below the Sommerfeld resonances and the perturbative annihilation formula!^ is adequate. 
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FIG. 2: Color factors for the different types of DM search experiments. The different multiplicity 
factors can be easily understood from the color flow in the figure. For direct detection the color of 
the incoming dark matter is the same as of the outgoing and so there is no multiplicity factor. For 
indirect detection it is an annihilation diagram, so just as for the thermal relic calculation there 
is a 1/2A^ suppression because a DM particle can only annihilate if it finds the anti-particle with 
the right anti-dark color. For colliders there is an 2N enhancement because any of the different N 
colors can be created and an extra 2 from Dirac vs Majorana. 

the mass range of interest (aronnd 1 TeV) it is above the nentrino backgronnd and potentially 
within reach of fntnre experiments. 


B. Indirect detection 

The annihilation cross-section relevant for indirect detection is snppressed by a l/2iV 
factor. Taking into acconnt Sommerfeld enhancement the cross-section is fnrther rednced 
relative to the standard “wino” model. This is becanse there is a resonant Sommerfeld en¬ 
hancement from weak interactions for dark matter masses in the 2 — 3 TeV range, bnt the 
enhancement is mnch smaller for masses aronnd 1 TeV. The “wino” model has been inves¬ 
tigated recently [3T1 - I35] . and is strongly disfavored by H.E.S.S. data. Onr SU{N)d model is 
not yet constrained by either H.E.S.S. or Fermi data for any N > 2, bnt the annihilation 
cross-section is close to H.E.S.S. sensitivity as shown in Fignre It is worth noting that 
the limits shown in Fignre [^assnme a specihc NEW prohle and there is a large nncertainty 





in these limits due to our limited knowledge of the dark matter distribution in the center of 
the galaxy. Also shown in the hgure is the projected reach of CTA [Sn], assuming 5 hours 
of observation time. One can see that CTA should have enough sensitivity to discover or 
rule out our model up to at least A^ = 10, and therefore is the most promising search for 
discovering dark matter in our model. 



FIG. 3: Indirect detection constraints from gamma ray line searches from H.E.S.S. (blue region) 
and projected sensitivity for gamma ray lines at CTA (gray region) assuming an NFW profile 
for the dark matter distribution in the galactic center, both taken from [33]. The red dots are 
the cross sections for dark matter annihilation into gamma rays in our models for A = 2 up to 
N = 10, which were obtained by appropriately rescaling the NLL cross section with Sommerfeld 
enhancement from |34) . For comparison we plotted in black the cross-section for the annihilation 
cross-section to photons in the standard “wino” model as a function of mass, also from |34j . 


C. LHC and future collider searches 

The multiplicity factor enhances sensitivity of collider searches to our DM in two ways. 
Most importantly, the predicted DM mass from thermal freeze-out is lowered, and thus 
a lower partonic center of mass energy is required to pair produce DM at a collider. For 
example, at the 14 TeV LHC and for DM masses near Mdm ~ 500 GeV the DM cross section 
scales as (l/M^j^)® because of the strong energy dependence of the parton luminosities. In 
addition, the cross-section for pair producing dark matter is enhanced by the dual state 
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multiplicity factor 2N. 



FIG. 4: Expected significance of missing energy (MET) searches for DM at the LHC and a future 
100 TeV collider. The solid black lines in each plot correspond to the sensitivity of the collider 
to “wino”-like DM, an SU{2)ujeak triplet Majorana fermion. The colored dots labeled by different 
A^-values correspond to our models in which the DM is a Dirac fermion with multiplicity N and 
mass chosen to yield the correct abundance from thermal freeze-out. 

In the left panel of Figure |^we show the expected sensitivity of the high luminosity LHC to 
the DM in our model. The solid dots correspond to DM with multiplicity N and Dirac masses 
chosen so that the correct thermal DM abundance is obtained. We see that the 14 TeV LHC 
is sensitive to DM with N > 7. For comparison, we also show the expected sensitivity to 
a standard “wino” SU{2)weak triplet with a Majorana mass. Note that the “wino” has 
the correct thermal abundance only for M^jino — 3 TeV. Existing mono jet searches from 
ATLAS [3Z] and CMS [3H] with 8 TeV collisions and a luminosity of 19.5 fb“^ already rule 
out N > 20. 

In the right panel we show that a 100 TeV future collider can discover our DM above 
backgrounds even for the smallest non-Abelian multiplet of V = 2 and perhaps may be able 
to rule out a Dirac “wino”. The signihcances for these plots were determined from parton 
level signal and background events which were computed with MadGraph [321 • The main 
irreducible backgrounds are due to jets plus Z or W with MET from decays to neutrinos. We 
computed signal and backgrounds to leading order (agaw) and assumed that the experiments 
will be able to limit background systematic uncertainties to 2%. For the “wino”, our results 
are consistent with the more sophisticated studies in ([iQlllI])- These references also showed 
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that a MET plus “disappearing track” search for the production of can improve sensitivity 
because it is free of irreducible SM backgrounds. 

Another observable consequence of this model is a change in the running of the EW 
gauge coupling WS- The multiplicity of DM leads to an 2N enhancement factor in the DM 
contribution to running of aw at one loop, which would be observable at the proposed 100 
TeV hadron colliders. 


IV. DARK GLUONS AS DARK RADIATION 


In this Section we turn our attention to the cosmic evolution of the energy density in dark 
gluons and its effects on the CMB. An important parameter which determines the effects of 
dark gluons on the CMB is the ratio of temperatures of the dark gluon plasma Td over the 
photon temperature T. 


Dark gluons are coupled to the thermal bath of SM particles through their couplings to 
the DM which is in equilibrium with the SM in the early universe. If ad is not too small 
the dark gluons also equilibrate with the SM plasma at early times, so that Td = T before 
DM freeze-out. To determine the smallest possible coupling ad for which the dark gluons 
are in chemical and kinetic equilibrium with the SM we consider the Feynman diagram in 
Figure At temperatures higher than M^, the thermally averaged cross-section for this 
process times the DM number density is given by 


n^{av) ~ 


TiawCtd 

'J'2 


(3) 


Comparing this to the Hubble rate H we hnd that the dark gluons are in equilibrium with 
the SM at T ~ for 

<^d> -TV-• 

aw Planck 

Thus even for ad as small as 10“^^ the dark gluons come to chemical and thermal equi¬ 
librium with the SM at temperatures of order M^. However, when the universe cools below 
T ~ dark matter becomes non-relativistic and its number density drops exponentially. 
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FIG. 5: The process through which dark gluons maintain equilibrium with the DM and the SM 
plasma for small ad- 

Then the rate for the process in Figure]^ becomes® 

( 4 ) 

and the dark gluons decouple at temperatures of about for ad 10 to about M^/20 for 
ad ~ 10“^. Below the decoupling temperature the dark radiation fluid evolves independently 
with a temperature Td which redshifts as 1/a. 

The temperature of the photon fluid also redshifts as 1/a for most of the universe’s 
evolution. However, when massive SM particles become non-relativistic they annihilate 
into the remaining lighter SM particles which effectively heats up the photons compared 
to the dark gluons (similarly to what happens to photons and neutrinos after neutrino 
decoupling m)- The ratio between the photon and dark gluon temperatures can be easily 
calculated in the instantaneous decoupling approximation by requiring that the entropy per 
co-moving volume is conserved independently in each fluid, 

ryd' 

^7 \9lJ 

where gl is the number of effective degrees of freedom in the SM plasma at the time of dark 
gluon decoupling and g( the number of effective degrees of freedom at any later time. 

The CMB places strong constraints on the energy density in relativistic particles at the 
time of recombination. This constraint is usually presented in terms of the number of 

® Note that for temperatures below the W mass, the W boson in Figure [^should be replaced by a photon 
and the aw in Eq. I by f^em- 
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effective neutrino species, Neff. The contribution of the dark gluons to Neff is given by 

(6) 

where the — 1 is the number of generators of SU{N)d and Tj, is the neutrino temperature. 
The ratio can be calculated using Eq. right at neutrino decoupling when neutrino and 
photon temperatures are still the same. Assuming that the decoupling between dark gluons 
and the SM happens at temperatures around 50 GeV one finds 

ANeff = 0.07 {N^ - 1). 

The strongest constrain on Neff comes from the 2015 Planck data [Tj, which found Neff = 
3.15 ± 0.46 at 95% conhdence to be compared with the SM prediction [33] N^j^ = 3.046. 
We see that this rules out > 4 and that = 3 is within the 2cr allowed range. However, 
we note that the Planck analysis assumes the ACDM model with one additional parameter, 
Neff. This limit could potentially be relaxed in our scenario where there are important 
differences as we will now discuss. 

The effects of dark radiation can be divided into so-called background effects and per¬ 
turbation effects. Background effects are due to a change in the average energy density 
in relativistic degrees of freedom and are not sensitive to any other properties of the dark 
radiation fluid. The largest background effect of extra radiation from relativistic degrees 
of freedom is to change the redshift of matter-radiation equality Zeq. Since Zeq is very well 
measured a fit to the data is forced to maintain the redshift of matter-radiation equality 
by simultaneously increasing the dark matter density (see e.g. [331 - I3U] ). This change in the 
matter density in turn requires a change in the Hubble parameter today Hq, in order to keep 
Qrn fixed (the ratio between the matter density and the critical density). Thus we see that 
a £t to the CMB data alone has an approximate flat direction in which an increase in A^e// 
can be compensated for by simultaneous increases in pdm and Hq. 

The perturbation effects are due to perturbations in the dark radiation fluid and thus 
sensitive to properties of dark radiation. In particular there are two additional parameters 
which distinguish different types of dark radiation, see e.g. [19]: the effective sound speed 
and the viscosity speed The dark gluons are a relativistic fluid and have =1/3, the 
same as for neutrinos. However—unlike neutrinos—the dark gluons have self-interactions 
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which come from the non-Abelian gauge kinetic terms. If the rate of dark gluon-gluon 
scattering is large compared to the Hubble parameter then the dark gluons are well described 
as an ideal fluid instead of a free-streaming fluid as for neutrinos. For an ideal fluid one 
has Cyjg = 0 instead of 1/3 as for neutrinos. The interaction rate between dark gluons is 
approximately given by 


— 1 ‘2rfi 

T ~ a^Td, 


(7) 


and one sees that as long as ad ^ 10“^^ this rate is larger than H during recombination. Thus 
on the time scale set by Hubble the dark gluons behave as a perfect fluid. Qualitatively, the 
interactions reduce the damping of overdensities from relativistic particles streaming out of 
gravitational potential wells. This has the effect that the CMB peaks are not as suppressed 
as they would be in the case of additional free-streaming dark radiation. 


The Planck Collaboration has performed a fit for Cg^ and with Neff fixed to the SM 
value 3.046, i.e. no additional dark radiation, and found that the parameters were in perfect 
agreement with the expected value for neutrinos, Cgg = c^jg = 1/3. However, Planck has not 
yet performed a fit for additional radiation ANeff which is allowed to have non-standard 
values for Ceg and Cvis- In particular, if future CMB experiments find evidence for a non-zero 
ANeff then measuring Cvis of this extra component would allow one to distinguish between 
dark gluons and free-streaming dark radiation, like dark photons or sterile neutrinos [men]. 


V. DARK MATTER-DARK GLUON INTERACTIONS AND LARGE SCALE 
STRUCTURE 

In this Section we study how interactions between dark gluons and DM affect the linear 
evolution of DM overdensities. The interactions have two important effects. One is the 
transfer of kinetic energy (i.e. temperature) from the DR to the DM. The other is that DM 
particles moving through the DR fluid experience a drag force. This drag slows the growth 
of large scale structures through gravitational clustering. Both rates can be computed by 
considering scattering of dark gluons with DM particles. The process is a generalization 
of Compton scattering to the non-Abelian case. The most important new feature is that 
scattering is dominated by the t-channel diagram shown in Figure which is divergent for 
small angle scattering and which only exists for non-Abelian gauge bosons. As we will see. 
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this interaction leads to interesting signatures in large scale structure which distinguish our 
scenario from other models with interactions between dark radiation (DR) and DM (see e.g. 
Dsmz]). In particular, for aa ^ 10 our model predicts a smooth suppression of the 
matter power spectrum at all scales which could resolve the conflict between the indirect 
measurement of as from Planck [1] and the direct measurement from BOSS [6]. 


P 


k 



P' 


k' 


FIG. 6: t-channel scattering of dark matter with dark gluons. 


To simplify the calculation of the energy transfer rate we consider the limit in which 
the DM temperature is negligible so that we can take the DM particles to be at rest. We 
then compute the rate of energy transfer to the DM from scattering [IHl EH] with a thermal 
bath of dark gluons. As in the well-known case of Coulomb scattering the cross section is 
dominated by small-angle forward scattering. To signihcantly impact the energy of a massive 
DM particle many collisions with the gluons are required. The collisions are uncorrelated 
so that the resulting momentum of the DM particle performs a random walk with 


E = 


p 


2M^ 2M^ 




2Mv 


Here 6p is the typical momentum transfer in a single collision and N is the number of col¬ 
lisions. After many such random scatters the resulting DM particle distribution is thermal, 
but not necessarily with the temperature of the gluon bath. The DM temperature depends 
on the relative size of the Hubble expansion rate and the energy transfer rate. 

The rate of energy transfer is calculated by averaging the energy transfer per collision 
over the initial Bose-Einstein distribution of the dark gluons 
d^k 


E = a 


(2vr) 


2{N^-l)f{k) X 


( 8 ) 


1 


d^k' 


d?p' 


,E,k J (2.r2k' . 

Here—and for the remainder of this Section—time derivatives are taken with respect to 
conformal time which is the origin of the scale factor a on the right hand side. Also, 
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f{k) = l/{exp{k/Tdr — 1) is the gluon thermal distribution function and it is multiplied by 
2{N‘^ — 1) for the spin and color of the initial gluon. The 1 + f{k') hnal state factor accounts 
for stimulated emission. Finally, the color and spin summed and averaged matrix element 
(keeping only the t-channel) is given by 


|M|2 


1 

2 


u) 


•> 


where the 1/2 is the disappointingly boring color factor of the t-channel diagram and s, 
t, and u are the Mandelstam variables. The integrals are straightforward to evaluate for 
p = 0) and give 


TT 


T3 

. -‘-di 


E = a{N -1) - log(a^ ) ^ = a- log(a^ ) 

O Mv TT 


Pdr 


(9) 


‘X '' TdrMx 

One subtlety one encounters is a logarithmic divergence at small t which stems from the long 
range of the interaction mediated by a massless gluon. In the thermal plasma of gluons the 
interaction range is made hnite by screening. This effect can be parametrized by including 
a Debye mass for the gluon, ~ d'd^dr ISH]- With the Debye mass the logarithmic 

divergence becomes the log(a^^) in Eq. This is analogous to the ubiquitous “Coulomb 
logarithm” in plasma physics. 


Using T = 2/3 E we obtain an equation for the evolution of dark matter temperature [51] 

p = -2^T,+a(/V^-l)^a^logK-‘)^(r*-T,) , (10) 

X 

where the Tdr — factor generalizes the energy transfer rate in Eq. to the case in which 
the thermal motion of the DM is not negligible. To understand the possible solutions to 
this equation one compares the size of the temperature transfer term with the Hubble red- 
shift term a/aT^ ~ aT-^T"^/Mpianck (during radiation domination). For large temperature 
transfer rates DM and DR are kept in thermal equilibrium by the interactions and = Tdr- 
Note that both terms scale with the cube of the temperature, thus if DM and DR are in 
equilibrium at one temperature they will be in equilibrium throughout radiation domination. 

For couplings ad < 10“® and assuming comparable starting temperatures ~ Tdr, 
the Hubble term initially dominates over the interaction term and the DM temperature 
plummets relative to the temperature of radiation: ~ versus Tdr ~ cl~^- This is 

the regime where our DM behaves like ordinary cold DM (CDM). However, note that the 
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Hubble term scales proportional to whereas the collision term scales like Tj^,. Thus for 

sufficiently small the two terms become comparable and the DM temperature switches to 
scaling proportional to a“^, keeping the ratio T^/T^r constant. In Figure]^ we show numerical 
solutions for the DM temperature as a function of scale factor for three representative values 
of ad- 



FIG. 7: Dark sector temperatures as a function of scale factor a. Shown are the dark gluon 
temperature (black dashed), and DM temperatures for three representative values of ad = 10“® 
(upper, blue), ad = 10“® (middle, red), and ad = 10“^*^ (lower, green). 

For values of ad for which DM and DR are in thermal equilibrium DM behaves very 
differently from ordinary CDM. The pressure from the dark gluons prevents the growth of 
DM overdensities during radiation domination. This can be seen as the sharp drop in the 
DM power spectrum in Figure For smaller values of ad the dark gluons still influence 
the evolution of DM overdensities, however the effects are more subtle and we employ the 
formalism of Ma and Bertschinger ISB to study them. 

Following [HI] we write down the linearized evolution equations for overdensities including 
the interactions between the DM and DR fluids in conformal Newtonian gauge. To avoid 
the complications of solving a full Boltzmann code we work with a simplihed scenario. 
We replace all relativistic energy density in the SM (i.e. neutrinos and photons) with an 
identical energy density which is made up of only photons. And we replace all matter (dark 
matter and baryons) with an equivalent energy density of only dark matter. In addition, 
we approximate by treating the photons as a perfect fluid (zero viscosity). This is only 
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true before recombination but soon after recombination photons contribute only a negligible 
amount to the energy density and thus to the evolution of DM overdensities. Since our goal 
is to demonstrate the effect of the interactions on the DM density perturbations we compare 
our scenario with interactions to the same scenario with ad = 0. 


In Fourier space the equations for the DM and DR overdensities are 


^DM — —doM + 3-^ 


OdM — - OdM + OR {(^DR — Gdm) + k t/j 

a 

4 

5 DR = —-^6 DR + 4-^ 

n j 2^DR I ;2/ , ^ PDM -1/n Q \ 

Odr — k —^- \- k Ip + —ar^ [Odm — (jdr) 


4 PDR 


( 11 ) 


5^ — — 0^7 4 ” 
o 

9^ = + k'^'ip 

k'^'ip + 3-('tp + -'ilj] = 


2Ml, 




where the dots represent derivatives with respect to conformal time, px is the average energy 
density of fluid X and 5x and 6x are related to the overdensity and velocity divergence in 
fluid X. We have also set the two metric perturbations pj and 0 equal because we are treating 
the photons and dark radiation as ideal fluids (no anisotropic stress) and did not include 
neutrinos which have sizable anisotropic stresses. The interaction between dark matter and 
dark radiation is encoded in the momentum transfer rate r“ m- It is dehned as the change 
in momentum = —ar“^ Py. which a DM particle with momentum P experiences due to 
friction as it is moving through the dark gluon fluid. Microscopically, the friction arises from 
collisions between DM particles and dark gluons and to compute it we evaluate 


p = a 


(Pk 


f{k)- 


(Pk' 


<Pp' 


, (2.)3. 

where now the initial DM momentum p is non-zero and we expand to hrst order in p/M^^. 
Employing the same approximations as for the energy transfer rate we obtain 








2 

dr 




(13) 


We integrate the equations for the overdensities from a = 10 when all modes of interest 
are well outside the horizon, until a = 10“^. We use initial conditions corresponding to 
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adiabatic perturbations: 


<^7 — ^DR — — —2'0 — C{k ), (14) 

where the initial perturbations C{k) ~ are determined by the physics of inflation. We 
dehne the DM power spectrum equal to the square of the perturbations, P{k) = at 

scale factor a = 10“^. To focus only on effects of the coupling between DM and DR we 
form a ratio where we divide the power spectrum with interactions turned on by the power 
spectrum with aa = 0. Note that since the equations are linear the initial values for the 
perturbations, C{k), drop out in the ratio. 



k (Mpc^b 

FIG. 8: Power spectrum including the DM-DR interactions normalized by the power spectrum 
with interactions turned off. The black dotted curve corresponds to = 10“®, the green dashed 
curve corresponds to ad = 10“®'® and the red line corresponds to ad = 10“®. The power spectra 
are defined proportional to at a = 10“®. The vertical yellow band labeled keg indicates modes 
which enter the horizon at matter-radiation equality, modes which enter the horizon earlier are to 
the right (larger k). The blue band labeled as indicates modes which the observable ag is most 
sensitive to. 

The ratios of power spectra for different values of ad are plotted in Figure]^ For the plot we 
chose the number of colors N = 2 and = 1.2 TeV. One sees that for ad = 10“® the power 
spectrum is strongly suppressed for modes which entered the horizon before matter-radiation 
equality. These are modes with k > keq ^ 0.015 Mpc“^. This should be expected because in 
this case the DM is in equilibrium with the DR bath throughout radiation domination. For 
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the smaller values of ad = 10“®'^ and 10“® the power spectrum is less affected with modes 
which entered the horizon earlier (larger k) suppressed more that those which entered later. 
Modes which enter the horizon after matter radiation equality are not suppressed for any of 
the couplings plotted. The light blue vertical band indicates the range of modes which the 
observable as is sensitive to (cxs is a measurement of the matter fluctuations in spheres of 
radius of 8h~^ Mpc). 

The smooth suppression of power at all scales that we are hnding is special to our scenario 
and stems from the fact that the momentum transfer rate scales with temperature as 
the same scaling as Hubble. Thus it is possible to arrange for the couplings to have a small 
effect but act over a large range of scales. This should be contrasted with cases where the 
interactions scale like a higher power of Tdr in which case they are important at high scales 
and have no effect at low scales. In such scenarios the power spectrum has a sharp cutoff 
at scales of the size of the horizon at the time when the interactions cease to be important, 
leaving larger scales unaffected and wiping out the smaller scales that entered the horizon 
at earlier times. 

The smooth suppression of power is exciting because it might help resolve two sources 
of tension in recent experimental results. A £t to the most recent Planck CMB data is 
used to £x the parameters of ACDM. Using the model, the Planck collaboration predicts 
agCDM _ 0 329 ± 0.014 (“TT + lowP”, Icr errors). This value is about 2a higher than 
“direct” measurements of large scale structure [U El [7] (for example, gravitational lensing of 
the CMB as measured by Planck gives = 0.802±0.012). A reduction of the predicted 

power spectrum due to DM-DR interactions as, it occurs in our model for ~ 10“®, removes 
this tension between the Planck £t and LSS data. Interestingly, the tension in cjg is currently 
also driving the Planck fit for N^ff and Hq to lower values, because in ACDM larger values 
for those parameters would correlate with even larger values for ug. Thus after including 
DM-DR interactions the Planck £t might prefer N^fj > 3 which would create more room 
for dark gluons. Furthermore, larger Neff is correlated in the ACDM fit with a larger value 
for Hq (to keep the position of the acoustic peaks in the CMB hxed [H IUHlG] ). This in turn 
would allow better agreement between the Hq values from Planck and supernova data [T], 
another area of mild tension in cosmological data. Clearly, a more quantitative analysis 
of this issue is desirable and requires including the non-Abelian dark radiation in a full 
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Boltzmann code. 


VI. CONCLUSIONS 


In this work we have studied a new type of dark sector with massless non-Abelian gauge 
bosons super-weakly coupled to the DM. There are many different possibilities for the DM 
coupling to the standard model. As an example, we chose our DM particle to transform as 
SU{2)weak triplets and fundamentals under the dark SU{N)d gauge group. 

Our model has three new parameters, the coupling constant aa, the mass of the DM 
particles M^, and the size of the gauge group N. Demanding the correct DM abundance 
from thermal freeze-out fixes the DM mass in terms of N, thus leaving a two dimensional 
parameter space. Constraints on this parameter space can be derived from several different 
experiments. 

The first set of constraints derives simply from the multiplicity of the dark matter and 
would even apply if the dark coupling constant were zero. The multiplicity of dark matter 
affects the usual WIMP searches (direct and indirect detection and collider searches). The 
effect is simply that dark color multiplicity factors enhance pair production and decrease pair 
annihilation. Therefore the DM mass required in order to predict the right DM abundance 
decreases by \/2N. It also increases the collider cross-section, placing this type of dark 
matter within easy reach of the proposed 100 TeV collider. The decrease in mass and 
in annihilation cross-section also removes the current tension between thermally produced 
SU{2)weak triplet dark matter and H.E.S.S data. Despite the decrease in annihilation cross- 
section, the SU{2)weak triplet DM model annihilation cross-section is within the projected 
reach for CTA. 

For ad ~ 10“^^ the dark gluons thermalize with the SM in the early universe and con¬ 
tribute to dark radiation. The limits placed on N^ff from Planck constrain N to be at most 
3. However, the self interactions of the dark gluons leave an imprint in the CMB which is 
distinct from that of free-streaming fluids like neutrinos or dark photons. This can be used 
to distinguish between the two types of radiation if future experiments establish the need 
for a non-standard contribution to 
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Finally we studied the effects of the interactions between DM and DR on the power spec¬ 
trum. We found that for > 10“® the interactions strongly suppress the power spectrum 
of modes entering the horizon before matter radiation equality and thus such couplings are 
ruled out. On the other hand, for ad ~ 10“®'^ the interactions predict a smooth decrease in 
the power spectrum, which can potentially solve the discrepancy between Planck and large 
scale structure data and the discrepancy between Planck and Supernova measurements of 
Hq. We hope to return to this in future work. 
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